Gravitational Waves versus Electromagnetic Emission in Gamma-Ray Bursts 
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The recent progress in the understanding the physical nature of neutron star equilibrium configu- 
rations and the first observational evidence of a genuinely Short Gamma-Ray Burst, GRB 090227B, 
allows to give an estimate of the gravitational waves versus electromagnetic emission in a Gamma- 
Ray Burst. 



We first recall that we have recently proved 1-3] how 
the consistent treatment of neutron star equilibrium con- 
figurations, taking into account the strong, weak, elec- 
tromagnetic, and gravitational interactions, implies the 
solution of the general relativistic Thomas-Fermi equa- 
tions, coupled with the Einstein-Maxwell system of equa- 
tions. This new set of equations supersed the tra- 
ditional Tolman-Oppenheimer-Volkoff (TOV) equations 
0, H1 1 which impose the condition of local charge neu- 
trality throughout the configuration (see e.g. Q and ref- 
erences therein). 

The solution of the Einstein-Maxwell-Thomas-Fermi 
coupled differential equations leads to a new structure 
of the star (see |l|, for details): the positively charged 
core at supranuclear densities, p > p nuc ~ 2.7 x 10 14 g 
cm -3 , is surrounded by an electron distribution of thick- 
ness > h/(m e c) and, at lower densities p < p n uc, a neu- 
tral ordinary crust. The equilibrium condition given by 
the constancy of the particle Klein potentials leads to 
a discontinuity in the density at the core-crust transi- 
tion and, correspondingly, an overcritical electric field 
- (m,/m e ) 2 B Cl where E c = m 2 e c 3 /(eh) ~ 1.3 x 10 16 
Volt/cm, develops in the boundary interface; see Fig. Q] 
In particular, the continuity of the electron Klein poten- 
tial leads to a decreasing of the electron chemical poten- 
tial p e and density at the core-crust boundary interface. 
They reach values p,f ust < n c ° vc and p C mst < Pcorc at 
the edge of the crust, where global charge neutrality is 
achieved (see Fig. [I}. We shall adopt some features of 
these neutron stars computed using the NL3 parameter- 
ization 0] of the phenomenological a-u>-p nuclear model 
[1]; we refer to [l[ for details. 

For each central density there exists an entire family 
of core-crust interface boundaries and, correspondingly, 
a family of crusts with different mass M crust and thick- 
ness Ai? crust . The larger p CTUS t, the smaller the thickness 
of the core-crust interface, the peak of the electric field, 
and the larger the M crust and Ai? crust . The configuration 
with p crus t = Pdrip ~ 4.3 x 10 11 g/cm 3 separates neutron 
stars with and without inner crust. All the above new 
features lead to crusts with masses and thickness smaller 
than the ones obtained from the traditional TOV treat- 
ment. The mass-radius relation obtained in this case 
have been compared and contrasted with the one ob- 
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FIG. 1: Upper panel: particle density profiles in the core-crust 
boundary interface, in units of cm" 3 . Middle panel: electric 
field in the core-crust transition layer, in units of the critical 
field E c . Lower panel: density profile inside a neutron star 
with central density p(0) ~ 5p nuc . We compare and contrast 
the structural differences between the solution obtained from 
the traditional TOV equations (locally neutral case) and the 
globally neutral solution presented in |l| . In this example the 
density at the edge of the crust is p C rust = Pdrip = 4.3 x 10 11 
g/cm 3 and \ a = h/{m a c) ~ 0.4 fm denotes the <r-meson 
Compton wavelength. 



tained from the locally neutral TOV approach; see Fig. [5] 
and [if for details. 

In Fig. [2] we show how our new neutron star the- 
ory is in agreement with the most up-to-date strin- 
gent observational constraints to the mass-radius rela- 
tion of neutron stars, that are provided by the largest 
mass, the largest radius, the highest rotational frequency, 
and the maximum surface gravity, observed from pulsars 
0. They are imposedby the mass of PSR J1614-2230 
M = 1.97 ± 0.04Mq_M, the lower limit to the radius 
of RX J1856-3754 [liTdotted-dashed curve), the 716 
Hz PSR J1748-2246ad (dashed curve), and the sur- 
face gravity of the neutron star in the Low Mass X-Ray 
Binary X7 from which 90% confidence level contours of 
constant Roo can be extracted [HI (dotted curves); see 
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1] for further details. 

The above constraints strongly favor stiff nuclear equa- 
tions of state such as the ones obtained from relativistic 
mean field models, which provide high maximum masses 
for neutron stars [9|. In addition, the radius of a canon- 
ical neutron star of mass M = 1.4M Q is highly con- 
strained to the range R > 12 km, ruling out a strange 
quark hypothesis for these objects. Our new neutron 
star mass-radius relation fully agrees with all the above 
requirements, for instance, we find that a canonical neu- 
tron star with M = 1.40M Q has a radius R = 12.31 km, 
for the NL3 parameterization of the nuclear EoS (see [l[ , 
for details). 
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TABLE I: Properties of GRB 090227B: Egf B is the total 
energy emitted in the GRB, B is the Baryon load, Ytr is the 
Lorentz factor at transparency, the cosmological redshift is 
denotes by z, the intrinsic duration of the GRB is At, and 
the average density of the CBM is {ticbm)- We refer to pi| 
for additional details. 
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FIG. 2: Constraints on the neutron star mass-radius relation 
(see [i| and references therein). We compare and contrast 
the theoretical M-R relation of globally neutral neutron stars 
(blue curve) obtained from the solution of the Einstein- 
Maxwell- Thomas-Fermi equations and locally neutral neutron 
stars (red curve) obtained by solving the TOV equations. Any 
mass-radius relation have a maximum mass larger than M = 
1.97±0.04Mq and should pass through the area delimited by 
the solid, dotted-dashed, dashed, and dotted curves. 

We now turn to the observations of GRB 090227B 
(see [lj| for details). The progress obtained from the 
Fermi-GBM and Konus-Wind satellites has been used to 
identify the new class of genuinely short GRBs: short 
bursts with the same inner engine of the long GRBs but 
endowed with a severely low value of the baryon load, 
B = M B c 2 /Ef f- B < 5 x 10~ 5 , where M B is the mass 
of the baryons engulfed by the expanding ultrarelativis- 
tic e + e~ plasma of energy E^f B . The emission from 
these GRBs mainly consists in a first emission, the peak 
GRB (P-GRB), followed by a softer emission squeezed 
on the first one. The typical separation between the two 
components is expected to be shorter than 1-10 ms. 

A special case is GRB 090227B. From the 16 ms time- 
binned light curves a significant thermal emission in the 
first 96 ms, which has been identified with the P-GRB, 
has been found 16]. The subsequent emission is identi- 



fied with the extended afterglow. The P-GRB of 090227B 
has the highest temperature ever observed, k B T = 517 
keV, where k B is the Boltzmann constant. The results of 
the fit of the light curve and spectrum of GRB 090227B 
are summarized in Table |TJ In particular we show the to- 
tal energy emitted EfJ^ B , Baryon load B, Lorentz factor 
at transparency r tr , cosmological redshift z, intrinsic du- 
ration of the GRB emission At, and average density of 
the CircumBurst Medium (CBM) (ticbm)', we refer to 
[l(| for further details. 

The above quantitative results lead to the conclusion 
that the progenitor of GRB 090227B is a neutron star 
binary: (1) the natal kicks velocities imparted to a neu- 
tron star binary at birth can be even larger than 200 kin 
s _1 [2~0 ] and therefore a binary system can runaway to 
the halo of its host galaxy (see e.g. fUj]), clearly point- 
ing to a ver y lo w average number density of the CBM 



(see e.g. [17|-|20|); (2) the very large total energy, which 



we can indeed infer in view of the absence of beaming, 
and the very short time scale of emission point again to 
a neutron star binary; (3) as we shall show below the 
very small value of the baryon load is strikingly consis- 
tent with two neutron stars having small crusts, in line 
with the recent neutron star theory [l[. 

We now infer the binary component parameters. It is 
clear that the merging of two neutron stars will lead to 
a GRB if the total mass of the binary satisfies 



Mi + M 2 > M c 



2.67M, 



o ) 



(1) 



where M cr it is the critical mass over which a neutron star 
undergoes gravitational collapse to a black hole. The 
numerical value reported in Eq. ([lj has been taken from 

a. 

Assuming for simplicity a binary with twin components 
Mi = M 2 = M, we obtain masses M = 1.335M© and 
correspondingly radii Ri = Rz = 12.24 km (see Fig. 2 
and [l|). The mass of the corresponding crust of each 
component is M crust ~ 3.6 x 1O _5 M0 and the thickness 
of the crust is Ai? crust ~ 0.47 km. 

The location of the binary in the very low interstellar 
density medium of galactic halos makes possible to probe 
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FIG. 3: Baryon load expected to be left by a binary neutron 
star merger, given by Eq. ([2}, as a function of the total mass 
M of globally (lower panel, units 10 -5 ) and locally neutral 
(upper panel, units 10 -2 ) neutron stars, for the case of GRB 
090227B. We have indicated the observed baryon load of GRB 
090227B, B — 4.13 x 10" 5 ; see Table Hand 



the neutron star theory and equation of state through the 
knowledge of the baryon load B inferred from the fitting 
of the GRB light curve and spectrum. The baryonic mat- 
ter which the GRB interact with is in these systems pro- 
vided by the material of the neutron star crusts ejected 
during the binary coalescence. Thus, a theoretical expec- 
tation of the baryon load B left in a binary neutron star 
merger is 



B = 



E 



e? bb 



where 



E 



2GM mTB M, 



corc Jw crust 



(2) 



(3) 



is the mass-energy of the crust material that will be po- 
tentially ejected; the factor 2 comes from the fact that 
both neutron star components contribute to the baryon 
load. 

In Fig. |3] we have plotted the theoretical baryon 
load given by Eq. © for GRB 090227B, namely using 
E? f B = 2.83 x 10 53 erg, as a function of the mass M of 
the globally and locally neutral neutron stars shown in 
Fig.H 

The agreement of the observed baryon load of GRB 



090227B (see Table U and [16|) with the low mass of 



the crust obtained from the globally neutral neutron 
stars of [l[ is evident (see Fig. [3]). It can be com- 
pared and contrasted with the ones obtained enforcing 
the local charge neutrality condition. For the specific 
binary neutron star system studied here we obtain a 
theoretical prediction of the baryon load from Eq. ©, 



B ~ E B ust /Egf B ~ 7.6 x KT 5 , or a mass of the 
baryons M B = E B ust /c 2 ~ 1.2 x lCr 5 Af , to be con- 
fronted with the one obtained from the fitting procedure 
of GRB 090227B, B ~ 4.13 x KT 5 , corresponding to 
M B = BE^ B /c 2 ~ 0.7 x Kr 5 M . The above theo- 
retical predictions of the neutron star crust mass M crust 
and consequently the value of E B ust and B have been in- 
ferred for a crust with a density at its edge equal to the 
neutron drip density pdrip ~ 4.3 x 10 11 g cm -3 . Neutron 
star crusts with densities/? < pdrip are predicted by the 
new neutron star theory [l| , there is still room for smaller 
values of the baryonic matter ejected in a binary process, 
and consequently to still shorter genuinely short GRBs. 

The mass-energy of the baryon ejecta obtained from 
the estimate ((2]) gives for locally neutral neutron stars 
values 10 2 -10 3 bigger than the ones analyzed before (see 
Fig. [3]), due to the more massive crusts obtained from 
the TOV-like treatment (see [l[ for details). It implies 
that Eq. © gives in such a case M B ~ 1CT 3 -1O" 2 M , 
in line with previous results obtained from the numeri- 
cal simulation of the dynamical evolution of neutron star 
binaries (see e.g. 0, 0|), where locally neutral neutron 
stars are employed. 

We are considering a neutron star binary with com- 
ponents (Mi,Ri) and (M 2 ,i?2) orbiting each other in a 
circular orbit of radius a (separation between the two 
neutron star centers), orbital angular velocity 



and total energy 



E 



G{M l + M 2 ) 



1 GM 1 M 2 

2 a 



(4) 



(5) 



The leading term driving the loss of energy via gravi- 
tational wave emission is given by 



dE _ 32 G 4 (Mi 
~dt ~~ 



M 2 )(M 1 M 2 ) 2 



5 c 5 



(6) 



which leads to a decreasing of the separation a with 
time and consequently a shortening of the orbital period 
P = 2it/lu dictated by P~ x dPjdt = {3/2)a,- 1 da/dt = 
-(3/2)E- 1 dE/dt Additional contributions to the 

gravitational wave power due to higher multipole mo- 
ments of the components such as angular momentum 
J and quadrupole moment Q (deformation) are con- 
ceptually relevant corrections to the above formula (see 
e.g. (2(| and references therein, for details) although 
quantitatively negligible. For instance, the first correc- 
tion due to the spin angular momentum J of the neutron 
star components is given by — 11/4 jujM in geometric 
units, where j — cJ/(GM 2 ) is the dimensionless angular 
momentum parameter. This correction is only of order 
10~ 2 for a binary orbit of very high angular frequency 
~ kHz and for neutron stars with M = 1.335M and 
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j = 0.4. We recall that the fastest observed pulsar, PSR 
J1748-2246ad, has a rotation frequency of 716 Hz [TBI) . 
which gives j ~ 0.51 {Mq/ M Q ) 2 = 0.26 J 45 with the 
latter value for a canonical NS of M = 1.4M Q , J45 is the 
moment of inertia in units of 10 45 g cm 2 . The first cor- 
rection due to the quadrupole deformation multipole mo- 
ment Q of the neutron star, given by -2Qw 4 / 3 M- 5 / 3 , is 
of order 10~ 3 for the same parameters with Q ~ 4 x 10 43 
g cm 2 ~ 3 km 3 , the latter value in geometric units. 

The gravitational waves emission dominates the energy 
loss during the spiraling phase while the electromagnetic 
radiation dominates from the coalescence with the final 
emission of a short GRB if the total mass of the binary 
exceeds the critical mass for neutron star gravitational 
collapse ~ 2.67M Q , see Eq. (JlJ. Thus, an upper limit 
for the gravitational wave emission radiated away can be 
obtained from the energy difference between the initial 
binary at time to = with separation ao and energy Eq, 
and the binary at time tf and separation a/ = R± + R%, 
with energy Ef, when the two components touch each 
other. An absolute upper limit for the gravitational wave 
energy emission, AEq^, can be therefore determined by 
the assumption of an infinite initial separation ao — > 00, 
i.e. 



a nmn 



E(t f ) - E(t ) 



1 GM\M% 

2 Ri + R 2 



(7) 



For the neutron star binary discussed in this work for 
GRB 090227B, we obtain the absolute upper bound 



1±E GW 



1 GM 



4 R 



9.6 x 10 52 erg = 0.054Af Q c 2 , (8) 



which in the case of the genuinely short GRB 090227B 
is one order of magnitude smaller than the emitted elec- 



tromagnetic energy E. 



GRB _ 



2.83 x 10 53 erg (see Table 



Uand 16]). It is also worth mentioning that indeed this 
numerical value (j5]) limits from above the results of full 
numerical integrations of the gravitational wave radia- 
tion emitted in the neutron star binaries during the entire 
process of spiraling and merging (see e.g. [23J). 

We have shown that the observations of the genuinely 
short GRB 090227B lead to crucial information on the 
binary neutron star progenitor. The data obtained from 
the electromagnetic spectrum allows to probe crucial as- 
pects of the correct theory of neutron stars and their 
equation of state. The baryon load parameter B ob- 
tained from the analysis of GRB 090227B, leads to a 
most remarkable agreement of the baryonic matter ex- 
pected to be ejected in a neutron star binary merger and 
validate the choice of the parameters of the binary com- 
ponents, M 1 = M 2 = 1.34M , and J?i = R 2 = 12.24 km. 
This represents a test of the actual neutron star param- 
eters described by the recent developed self-consistent 
theory of neutron stars [H that takes into account the 



strong, weak, electromagnetic and gravitational interac- 
tions within general relativity. In view of the above limits 
given in Eq. jSJ), it is clear that the emission of electro- 
magnetic radiation in a GRB by a binary neutron star 
system is at least one order of magnitude larger than the 
gravitational wave emission. 
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